!!
module force
use mod_2dflu
implicit none
integer,dimension(2),save :: kforce
double precision :: famp
public
namelist /force_pars/ &
  famp,kf
contains
contains
!------------------------------
subroutine add_cont_force()
  integer :: i1,i2,ireal,iimag
  i2=1
  i1=kforce(1)+1
  ireal = 2*i1-1
  iimag=2*i1
  omega(ireal,i2)=omega(ireal,i2)+famp*kf/2.0d0
endsubroutine add_cont_force
!------------------------------
subroutine gen_force
  integer:: i1,i2
  real*8::x,y
!
  fomega = 0.0d0 ! Initialize the force field
!
! Generate the forcing field in real space
!
  do i2 = 1,n2
     y = i2*dy
     do i1 = 1,n1
        x = i1*dx
	fomega(i1,i2) = -famp*kf*dcos(kf*x)
     enddo
  enddo
!
! Fourier tranform the force field
!
!  call rfftwnd_f77_one_real_to_complex(pfor,fomega,0)
!
end subroutine gen_force
!---------------------------
end module force
